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Abstract 


This  paper  introduces  and  analyzes  two  ways  of  extracting  the  hydrostatic  pressure 
when  solving  Stokes  problem  using  the  p  version  of  the  finite  element  method.  When 
one  uses  a  local  Hl  projection,  we  show  that  optimal  rates  of  convergence  for  the 
pressure  approximation  is  achieved.  When  the  pressure  is  not  in  B$y  or  the  value  of 
the  pressure  is  only  needed  at  a  few  points,  one  may  extract  the  pressure  pointwise 
using  e.g.  a  single  layer  potential  recovery.  Negative  norm  and  interior  estimates  for 
the  Stokes  velocity  are  derived  within  the  framework  of  the  p  version  of  the  F.E.M. 
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1.  Introduction. 

There  are  cases  in  which  a  continuum  is  subjected  to  an  incompressibility  constraint, 
such  as  a  divergence  constraint.  We  will  investigate  some  problems  dealing  with 
stability  of  a  class  of  numerical  discretizations  of  such  problems.  As  an  important 
example  we  consider  Stokes  equations  in  two  dimensions: 

-AU  +  VP  =  F  in  fl  C  R2  ,  . 

V-U  =  0  in  n 

with  appropriate  boundary  conditions  on  dCt.  In  the  standard  weak  formulation, 

FindU€^C[ffl(n)f  and  P  €  W  C  L2(H)  such  that 

o(U,v)  +  6(v,P)  =  (F,  v)  VvGF  (1.2) 

ft(TJ,  q)  =  0  Vg€W 

where  the  bilinear  forms  a  and  b  are  given  by 

«(U,v)  =  2/nEi.,c,J(U)cl,(v)dx  /13x 

6(v,  P)  =  -  /n  V  •  vPdx  (  ' 

and  (F,v)  denotes  the  usual  [L2(fl)]2  inner  product.  e,,(v)  is  the  symmetric  part 
of  the  deformation  gradient  \  ( -f  ~<).  The  pair  of  spaces  (V,  W)  depends  on  the 
boundary  conditons.  For  no-slip  boundary  conditions, 

{V,fV)  =  (|^(n))*,{«e  L*((l)  :  I  qdx  =  0}) 

■'n 

and  for  stress-free  boundary  conditions, 

(V,W)  =  ({v  €  H\ 0)  :  £lJ(v)  =  0>\L-(n)). 

We  discretize  by  choosing  finite  dimensional  spaces  VN  C  V,WN  C  W ,  such  that  we 
are  to 


Find  Vjy  and 
a(U/v,  v)  +  6(v,  Pat) 

&(Un  ,q)  = 


Pn  6  W/v  such  that 
(F,v)  Vv  €  VN 
0  Vg  6  WN 


(1.4) 


The  main  problem  is  whether  it  is  possible  to  select  a  sequence  of  pairs  (V^W^) 
that  is  stable  and  that  has  good  approximation  properties.  The  BabuSka  -  Brezzi 
condition 

min  max 
iewN\{o)vevN\{o) 


fn  V  •  v  q  dx 


>  c  >  0 


(1.5) 
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with  c  independent  of  N ,  is  sufficient  to  guarantee  stability  in  the  sense  that  the 
discretization  errors  are  within  a  uniform  constant  of  the  distances  from  (U,  P )  to 
(Vjv.W'n)  (cf.  [4J,[8]). 

We  shall  concern  ourselves  with  methods  using  high  degree  polynomials  to  ap¬ 
proximate  the  solution  to  the  Stokes  equations.  Among  possible  choices:  spectral 
methods  (cf.  [ll])  and  the  p-version  of  the  finite  element  method,  we  shall  deal  with 
the  latter.  Accuracy  is  improved  by  increasing  the  polynomial  degree  p  in  the  (usually 
small  number  of)  elements  in  the  mesh,  see  [5]. 

We  saw  in  [16]  that  the  pair  ([Qp]2,  V  ■  Qp)  was  unstable,  since  the  right  inverse 
(V-)"1  :  WN  — >  VN  had  a  norm  in  B{L2,Hl)  that  was  growing  linearly  with  p.  When 
WN  =  V  •  V/v,  it  is  well  known  that  at  uniform  bound  on  (V-)_1  in  B(L2,Hl)  is 
equivalent  to  the  BabuSka  -  Brezzi  condition  (cf.  [23]  ).  Similar  results  hold  for  other 
pairs  in  which  the  degrees  of  the  polynomials  in  the  (V)y,Ww)  pair  is  one  less  in  WN 
than  in  VN,  see  [16].  That  the  inf-sup  constant  tends  to  0  like  ^  as  p  — *  oo  does 
not  mean  automatic  suboptimal  covergence  for  the  velocity  as  was  explained  in  [16], 
e.g.  But  the  pressure  convergence  rate  can  be  slowed  down  by  1  as  was  seen  in  an 
example  there,  and  following  another  example  does  not  necessarily  have  to,  see  [16]. 
We  will  try  to  recover  the  pressure  through  some  additional  effort  -  solving  a  Poisson 
problem  -  at  optimal  rate. 

From  (1.1)  we  see  that  P  satisfies 

A  P  =  VF  in  0 
§f  =  (F-t-AU)-n  on  30 

and  fn  Pdx  =  0.  n  denotes  the  unit  outward  normal  defined  almost  everywhere.  (1.6) 
is  a  Poisson  problem  for  which  the  accuracy  of  any  numerical  procedure  will  depend 
on  the  accuracy  of  approximation  of  U. 

2.  Stability  and  velocity  convergence. 

A  lack  of  stability  (deterioration  of  the  inf-sup  constant)  does  not  imply  suboptimal 
convergence  of  the  velocities.  Let 

Zn  —  {v  €  Vff  :  v  =  0  on  3H,  V  •  v  =  0  in  0}  (2.1) 

where  we  have  taken  no-slip  boundary  conditions.  Let  3  V  •  V^.  Then  U/v  is  an 
elliptic  projection  of  U  onto  Z ^  such  that 

II  U  -  Ujy  ||i,n<  min  ||  U  -  v  ||i,n  (2.2) 
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where  |j  •  ||*,n=||  *  ||ff'(n)  the  standard  Sobolev  norm,  see  [l].  Therefore,  there  exist 
stream  functions  <p,<Pn  G  Hq(Q)  such  that 

U  =  V  x  <p,  UN  =  V  x  <pN  (2.3) 

where  V  x  <p  =  (- §J,  §f).  Let  R  =  (-1,  l)2  and 

Qp  =  span{xmt/n  :  0  <  m,n  <  p}  (2.4) 

be  the  spaces  of  polynomials  of  separate  degree  less  than  or  equal  to  p.  We  now  have 
for  one  square  element. 

Proposition  2.1  If  Vn  =  [ Qp  n  Hq(R)}2  and  W ^  =  V  •  VN,  then  for  p  >  4, 

U tv  is  an  elliptic  projection  onto  with 

II  U  -  U*  ||i,*<  Cp~M  ||  U  Hm+1.1?  (2-5) 

Proof:  Let  YN  =  ((Vx)_1(Zjy))  DHq(R).  Then 


min 

vezN 


U- v 


<  ||  V  x  {<p-  0)  ||lifl 

<  min  ||  <p  -  0  II 
1>€Vn 

<  Cp- M  ||  p  IIm+i.r 

<  Cp~M  ||  U  ||M+ipfi 


which  combined  with  (2.2)  and  the  observation  that  for  p  >  4,  Yn  =  Qp  n  HI  = 

(1  —  x2)2(l  —  y2)2Qp~4  0  and  has  optimal  approximation  properties,  using  an  inter¬ 

polation  argument  due  to  BabuSka:  Optimal  rates  of  approximation  were  established 
in  [25]  for  M  >  |  and  are  clear  at  M  =  0,  upon  which  it  follows  through  interpolation 
by  the  /C-method  for  intermediate  values.  The  H2  trace  constraints  interpolate  as 
expected,  as  shown  in  [15]. 

□ 

Remark  2.2  When  F  is  0  or  sufficiently  smooth,  the  rate  appearing  in  (2.5)  is  not 
optimal  for  typical  corner  singularities  that  appear  in  the  solution  of  elliptic  boundary 
value  problems  on  polygonal  domains.  It  follows  from  the  analysis  in  [6],  that  if  U  is 
of  the  form  ra|  logr|4^(0)  -  expressed  in  polar  coordinates  at  the  corner  -  and  <f>  is 
smooth,  then 

min  ||  U  -  v  ||,< 


which  is  near  twice  the  rate  of  (2.5)  as  U  G  H*  for  s  <  a  +  1.  Let  us  briefly  recap 
approximation  theoretical  results  for  the  stream  function  <p.  Let  <p  =  r7|  log  r|?<£(0) 
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(-7  =  a  + 1).  Combining  the  results  in  [12],  Chapter  7,  and  [9],  one  sees  that  for  corner 
angles  u  ~  .8128tt,  the  first  non  Gaussian  integer  pole  of  a  certain  transcendental 
equation,  which  determines  the  strength  7,  is  double  and  q  =  1.  In  [25]  it  was  shown 
that  then  the  following  approximation  result  holds 

min  ||  <p-rp  |]2<  C|  logpjp-2^-1).  (2.6) 


with  C  independent  of  p.  By  the  above  interpolation  argument  one  can  get  |  log  p\e 
for  any  0  <  0  <  1  and  C  independent  of  6.  This  indicates  that  the  log  factor  is  not 
important  in  (2.6).  The  singular  form  of  the  solution 

ip  =  r'7!  logr|?<£(0)  (2.7) 


can  also  stem  from  the  specification  of  Dirichlet  data  for  U  on  the  boundary.  We  shall 
mention  an  example  of  this  later  and  note  now  that  in  such  a  case,  the  |  logp|  factor  is 
not  present  iff  it  is  not  present  in  the  <p  satisfying  the  boundary  conditions  and  this  sin¬ 
gularity  is  stronger  than  that  inherited  from  the  geometry. 

In  view  of  the  previous  remark,  let  us  introduce  a  rate  of  approximation  function 
92  which  expresses  the  hitherto  known  results,  cf.  thm.  6.1  in  [25] 


mm 


V-i>  ll><  C!R(p,t,p) 


(2.8) 


where  p  is  a  regularity  index 
(0,0,  s) 


P{<P)  = 


(T»0il  +  Tf  -  *) 


so  that  we  may  define 


9? 


=  J  P(2_t) 

1  |  logpY 


V  ll< 


when  <p  £  H‘  \  H'+t ,  Ve  >  0 
is  assumed  known  only, 
when  p  £  \  Hi+ri,  Ve  >  0  and 

has  known  dominating  singularity: 
ip  =  r1!  logr[,</»(0)  in  local  polar  coordinates 
at  the  boundary  and  V  —  IZ'P  £  H* ,  s  >  27 


when  p(y?)  =  (0,0,  s),  t  <  s 


SPl’P  !(l  l)R(ip)  when  p(’p)  =  (-»,«,  1  +  1  -  e),  e  >  0 


(2.9) 


(2.10) 


R{ip)  involves  \\<p  -  an<*  the  li  norm  of  the  vector  of  coefficients  of  ip. 

Statements  similar  to  those  in  prop.  2.1  can  be  made  for  spaces  of  polynomials 
of  total  degree  <  p  and  those  used  in  PROBE  [28]  for  p  >  7.  These  are  defined  as 
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follows 

=  span{  xmyn  :  0  <  m,n  <  p} 

=  span{  x myn  :  0  <  m,n;  m  +  n  <  p  (2.11) 

or  (m  =  p  A  n  =  1)  or  (m  =  1  A  n  =  p)} 

=  Pp  ®  {xpy,xyp}, 

respectively.  The  latter  is  convenient  when  approximating  Dirichlet  boundary  condi¬ 
tions.  If  one  relaxes  the  no-slip  boundary  conditions  along  part  of  the  boundary,  the 
threshold  values  decrease  accordingly. 

([Po  ]2,  Pp~l)  -  or  any  other  pair  with  a  similar  degree  selection  -  is  a  natural  choice 
in  the  sense  that  it  reflects  the  approximability  of  the  velocity  -  pressure  pair  given  by 
their  regularity  and  that  it  was  proven  h-stable  for  p  >  4  when  excluding  exceptional 
meshes,  see  [23].  The  stability  does  not  seem  to  be  restored  by  merely  restricting  the 
pressures  to  lie  in  Pp~2.  [22]  reports  the  behavior  of  the  inf-sup  constant  as  p-1/*  for 
large  p.  In  fig.  1  we  give  the  smallest  positive  singular  value  of  the  discrete  Stokes 
matrix  for  ([Pq]2 ,  Pp~3) .  Again  p~ l!A  seems  to  be  the  asymptotic  behaviour. 

It  is  possible  to  stabilize  at  a  cost.  In  [13]  and  [14]  ,  etc.,  it  has  been  proposed 
to  add  weighted  least  squares  terms  to  the  energy.  The  problem  possesses  coercivity 
and  is  only  in  the  very  limit  of  saddle  point  nature.  However,  it  seems  that  the 
least  squares  terms  added  require  that  U  G  \H2\2  ,  which  is  more  regularity  than  we 
would  like  to  require  a  priori.  In  fact,  the  examples  computed  in  [16]  indicate  that 
H2  regularity  of  U  was  sufficient  for  the  pressure  to  be  approximated  optimally  by 
the  mixed  p  F.E.M.  We  will  give  a  brief  explanation  of  this  later.  One  way  to  get  a 
stable  (Vjv,VVyv)  pair  is  derived  in  the  following  example. 

Example  2.3  A  natural  idea  for  proving  a  uniform  lower  bound  on  the  discrete 
inf-sup  constant  is  to  attempt  to  inherit  the  continuous  bound  by  constructing  the 
maximal  right  inverse  of  the  Div  operator.  This  general  avenue  was  thoroughly 
explored  in  [30],  [23]  and  references  listed  there.  In  the  context  of  the  p  version  we  refer 
to  [16].  We  will  mimic  the  construction  in  the  nondiscrete  case  where  uniform  bounds 
exist  also  for  nonsmooth  polygonal  domains  under  Dirichlet  boundary  conditions,  [2]. 

Let  R  =  (-1,  l)2,  $N  C  Pp+1  D  #<}(£)  and  WN  =  A$N.  Clearly 

W* 


Pp 

Sp 


V- 
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with 

II  (V-)-1?  Hi  =  ||  (V.^AA-1?  II, 

<  II  VIA'1?)  ||, 

<  II  A-1?  ||2 

<  C  ||  ?  ||o  V?  6  L2(R), 

where  C  depends  only  on  R ,  not  on  q ,  giving  a  uniform  bound  on  the  operator  norm 
of  Div-1.  We  may  therefore  select 

WN  =  A(P0P+1)  and 

v„  =  V(P'+1) 

with  no  regard  to  boundary  conditions  and  under  no-slip  boundary  conditions,  we 
instead  set 

WN  =  A((N}0)]2PP-7)  and 

Vs  =  V([JVl0)]JPp-7) 

where  iv{0^  =  (1  —  x2)(l  -  y2)  is  the  first  internal  shape  function.  As  Div  1  was 
uniformly  bounded  in  B{L2,Hl),  the  pair  (V$^,A$//)  is  p-stable,  where  $N  = 

[Ar(°)]2pp-7_ 

□ 

Both  choices  seem  unnatural,  however,  from  the  regularity  point  of  view  as  well 
as  from  the  stand-point  that  you  would  like  to  be  able  to  combine  the  p  and  the  h 
versions. 

3.  Negative  norm  estimates 

For  future  purposes  we  shall  derive  error  estimates  in  negative  norms  for  the  velocities. 
As  seen  in  section  2,  it  is  possible  to  characterize  the  velocity  convergence  by  that  of 
stream  functions.  U/y  is  an  elliptic  projection  of  U  onto  Zn  so  that 

II  U  -  U*  ||/r-(n)=||  V  x  (<p  -  <pN)  ||«*(n)<||  <P~  <Ps  ||/f-+‘(n) 

with 

a(V  x  [<p  —  <Pn)>V  x  x)  =  0,  vXGy/v  =  p0a(n)n(Vx)-1zw 
c(<p-ipNtx)  =  0 
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where 


u(u,t/)  =  J 


n 

(if  stress-free  b.c.,  get  In(uxxvxx  +  uvvvuy  +  2 uxv 
bound  ||  e  ||_*  for  some  nonnegative  inteqer  k. 


AuAvdx 


Vu,  v  e  h2{ci) 

vxy)dx).  Let  e  =  ip  —  <pN.  We  need  to 


Proposition  3.1  Let  Cl  —  (  — l,l)2,fe  G  N  and  e  be  the  difference  between  the  exact 
((p)  and  the  discrete  (pn)  stream  functions.  Then 

II  e  ||-fc,n<  Cdt(p,t,p(<p))p-k-2  (3.1) 

where  9?  was  defined  in  (2.10). 


Proof: 


n  (e,t/) 

e  ||-fc,n=  sup  jp-Tr— 

«€C«(n)\{o}  ||  v  ||*>n 


Now  we  define  xp  to  be  the  solution  of 


A2ip  =  v  in  H 
0  =  =  0  on  dfl 


(3.2) 


v  ^  Co°(^)  =>  0  ^  C°°(n)  and  ||  xp  ||s,n<  C  ||  v  ||,-4,n  as  we  will  discuss  in  the 
remarks  following  the  proof. 


(e,  v)  =  (e,  A 20)  =  c(e,  xp)  =  c(e,  xp  -  i j)  Vp  G  YN 


such  that 


|(e,v)|  <  II  «  Ikn  min  ||  xp  -  p  |(2.n 

<  II  «  lU.n  Cp -*-1 1|  i,  |U+4,n 

<  Cp-*-1  II  e  ||,.„  II  V  |U,n 

<  Cp-‘-!»(p,(,p(p))||„||tin 


using  p(xp)  =  (0,0,4  +  k)  and  the  regularity  assumption  on  <p. 


□ 

Remark  3.2  It  is  not  trivial  that  (3.2)  can  be  solved  with  the  indicated,  regular  shift 
theorem  holding.  It  is  most  easily  seen  from  a  construction  found  in  Serbin,  [24|  or  [7], 
and  utilizes  the  geometry  in  an  essential  way.  Let  g,  denote  the  reflection  about  T,, 
the  four  straight  line  segments  forming  dCI,  i  =  1, ■  •  •  ,4.  Let  H,  =  g,n,t  =  1, •  •  •  ,4. 
Now  let  0  |n.  (x)  =  -rp{gx).  Since  gx  G  Cl  and  xp  G  C°°{Cl),xi  G  C°°(nuU'=1  H,).  Let 
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H  =  (— 3,3)2  and  define  xp  as  the  odd  extension,  once  more,  into  the  4  corner  regions 
(±1,±3)2.  Now  xp  G  C°°(fj)  and  since  A2  is  even,  xp  satisfies  (3.2)  on  0  with  v  the 
function  defined  through  the  same  odd  extension  procedure.  Now  standard  interior 
estimates  apply,  fl  CC  0  and  ,  easily,  3 C,  ||  v  j|s-4,n  =  C  ||  v  ||a_4  n !  II  ^  l|»,n<||  II,, n- 

□ 

Remark  3.3  The  above  remark  and  thus  proposition  hold  true  for  a  rightangled 
triangle:  fl  =  {(z,  y)  G  fi!2  :  0  <  a:  <  1  A  |  y  |  <  1 } .  Now  we  use  that  A2  is  not  only 
even  but  also  symmetric  (in  i  and  y).  (When  extending  xp  as  an  odd  function  about 
the  cathetae,  xpTX{x)  =  xpyv(gx),  e.g.).  The  same  conclusion  applies  to  the  Laplacian, 
as  is  well  known  for  the  square. 

□ 

Remark  3.4  In  general,  such  a  shift  theorem  will  hold  when  all  angles  are  sufficiently 
small,  [15].  In  [12]  Chapter  7,  it  is  shown  for  A2  that  each  corner  of  interior  angle 
a:  supports  singular  functions  with  strengths  commensurate  with  the  poles  t  of  the 
following  transcendental  characteristic  equation 

sinh2(ra)  =  r2  sin2  a. 

See  also  [18]  and  [20].  The  roots  of  this  equation  were  tabulated  in  [9],  It  has  been 
shown  that  for  a  G  (0, 7r),  the  first  root  has  imaginary  value  (for  this  particular  setting 
the  factor  determining  the  strength  of  the  singularity) 

dim  t  >  ~ 
a 

Therefore  the  first  singular  functions  belongs  to  Hs  when 

n 


We  list  the  bounds  for  two  cases  of  special  interest.  For  k  =  0  in  (3.1),  s  =  4  and  it 
is  sufficient  that  a  G  (0, 7r/2] .  For  k  =  1  in  (3.1),  s  =  5  and  a  G  (0, 7r/3]  suffices. 

□ 

Remark  3.5  Concludingly,  the  shift  theorem  for  (3.2)  and  hence  the  estimate  (3.1) 
will  hold  for  domains  with  smooth  boundary,  rectangles,  triangles  with  two  sides  equal 
and,  with  a  given  ceiling  on  the  shift’s  upper  Sobolev  index,  for  polygonal  domains 
with  sufficiently  small  angles. 

Corollary  3.6  If  <p  is  singular  on  d fl  xuith  q  =  0  in  (8.7),  i.e.  p[<p)  =  (7,0, 1  +  7  - 
e),  e  >  0,  and  k  G  N, 

II  <  IUn< 


(3  3) 
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4.  Interior  estimates 

We  will  first  derive  such  estimates  for  unions  of  elements  within  Cl.  By  a  union  of 
elements  we  mean  a  single  element  or  a  patch  of  elements.  Let  K0  CC  Ki  be  unions 
of  elements  of  the  triangulation  or  lattice  such  that  Ki  is  either  a  rectangle,  a  triangle 
with  two  sides  of  equal  length,  or  a  polygon  with  sufficiently  small  angles.  We  proceed 
using  the  duality  ideas  of  [19].  We  shall  for  simplicity  take  a  most  typical  case,  F  =  0. 
Also,  we  will  assume  p(<p)  =  ("7,0, 1  +  7  —  e),  Ve  >  0.  Obvious  modifications  can  be 
performed  for  other  cases. 

Lemma  4.1  Let  s  >  0  be  an  integer.  Then 

II  «  II-..K.S  C(p- <*+1>  ||  ,  H,.*,  +  ||  «  II...,,,.)  (4.1) 

where  C  depends  on  s,Ko,Ki. 


Proof:  Let  B  be  a  ball  separating  K0  from  K 1,  i.e.  K0  CC  B  CC  K 1  and  let 
u>  £  C$°(B)  with  u>  =  1  on  K0  (B  need  only  have  C°°  boundary).  Then,  for  s  >  0, 
we  have 

II  e  |U *,<||  we  11-1,^=  sup 

/€ffn*(K,)\{0}  II  7  IU.K, 

But  V/  €  H’iKi),  3lv  £  H‘+4(Ki)r\Hl{K\)  such  that  Cl(tj,v)  =  (rj,f),Vrj  £  H^{KX) 
with  ||  V  ||.+4l|f1<  C  ||  /  H,.*,.  Here  c^-,-)  denotes  c|*,(-,-).  Thus 

II  e  || C  sup  77 — jj - 


c1(we,v )  =  J  A(we)A vdx  =  Ci(e,wv)  -f  R(e,u,v) 


where 


R(e,w,v)  |  =  |  J  (eAwAu  +  2Vw  •  VeAu  -  AeAwv  -  2AeVw  •  Vv)dx  | 

=  1/  e|A  wA«  —  2V  •  (AvVw)  -  A(uAw)  -  2A(Vcu  •  Vi/)]<iz 
<  C(u>)  ||  e  ||_,_ji*1||  v  ||, +4.*, 
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using  Green’s  formulae.  Also,  we  use  the  fact  that  for  some  0  €  Yn(K  i), 

|ci(e,wv)|  =  \  Cl{e,uv  -  tf))  \ 

<  C  ||  e  ||2,*,||  luv  -  rp  Ha./r, 

<  C(w)/>-<*+1)  ||  e  lU./r,  II  v  ||.+,, /r, 

using  that  K\  is  a  union  of  elements.  Otherwise  we  would  get  ||  e  ||2,n-  This  concludes 
the  estimate  (4.1). 

□ 

Lemma  4.2  Let  (p  G  H2(Q),pN  G  Y N  with  p  >  2  and  let  k  G  IV.  Then 

II  e  ||o,/r„<  C(p~ 2  ||  e  Ik*,  +  ||  e  H-*,*,)  (4.2) 

tu/iere  C  depends  on  K0,Ki,k. 

Proof:  For  j  G  Z+,  denote  the  index  i(j)  =  Let  K0  CC  Ki/2  CC  •  •  •  CC  K^k)  = 
Ki  be  unions  of  elements.  From  the  previous  lemma  with  s  =  0,  we  get 

II  «  ||o,Ko—  C'(P-2  II  e  \U,k1/2  +  ||  e  ||_i,i<-1/a) 

Successively  applying  the  previous  lemma  to  estimate  ||  e  ||_;,jrl(y) ,  we  get  the  estimate. 

□ 

An  easy  consequence  of  (4.2)  is  obtained  by  interpolation  between  (4.2)  and  the 
obvious  ||e||2,*„  <  C{\\e\\2tKl  +  ||e||_fc(/fl)  to  get 

IMIi.tfo  <  c'(P_1|lelk«’i  +  ||cll-*.«-,)  (4-2') 

The  next  lemma  provides  a  local  growth  estimate  on  the  if2-norm  of  <Pn. 

Lemma  4.3  Let  k  G  N.  Then 

II  <Pn  lk*»<  C(p-min<1^-l)>  ||  <pN  ll2.iT,  +  ||  <Pn  H-fc./r,)  (4.3) 

where  C  depends  on  K0,Ki,k  and  <p. 

Proof:  Let  n  be  the  elliptic  projection  of  Hq{Ki)  onto  Y^(Ki)  defined  by: 


cj(v  -  xv, ip)  =  0,  V0  G  Yfj{K\). 
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n  satisfies 


™  H2.IT,  <  C\Jci(nvtnv) 


ci{nv,ip) 
sup  77-7-77 - 

^ei7v(i)i)\{0}  ||  ^  H2./T, 

ci(v,rp) 
SUP  77-77 

V>eKAr(/Tx)\<0}  ||  0  lk#r, 

<  c  II  U  II 2,/C  1  • 


<  c 

<  c 


Let  B  be  a  ball  separating  K0  from  K\,  i.e.  K0  CC  B  CC  KX)u  G  Cq3(B)  with  w  =  1 
on  Kq.  Then 


||  <PN  ||2,K„<||  U<Pn  Ik.iC,  <  II  UPN  -  '*{u(Pn)  ll2.iT,  +  II  ^{^Pn)  H2.iT, 

Let  p{p)  =  (7,0, 1  +  -7  —  e)  so  p  e  fl),  for  7  >  1.  Then  up  6  i3rl+',_£(fl)  and 

is  zero  in  a  neighbourhood  of  dKu  so  there  exists  ip  6  Y^(K j)  so  that 


\u<p-ip  ||2<  C(u,p)p  2(7  1}  ||  p  ||2 


Also  ||  p-  <pN  ||2<  C(<p)p  2<7  l)  ||  p  ||2.  Finally,  ||  p  ||2<||  p  -  pN  ||2  +  ||  <Pn  ||z< 
C{p)  ||  <Pn  || 2,  so  that 


||  upN  -  n(upN)  ||2iKl  <  C  ||  upN  -  ip  || 2,ir , 

<  C[\\u(p  -  PN)\\i,Ki  +  ||w^  -  V'lU.iT,] 

<  C(w,y>)p_2(7_1)  ||  <pN  ll2.iT,  • 


The  other  addend 


II  n{u<pN)  ||2,k,<  Cci(n(u><pN),<p)  =  Cci(uj<pN,0) 

where  <p  =  n(upN)/  ||  n(upN)  Ik*,  so  that  ||  ip  ||2.k,=  1. 

ci{u<Pn,p)  =  Ci(<pN,ojip)  +  R{pNfw,p) 

as  before  and  we  get  that 

I  Ct{upN,ip)  |<  C(p~J(7_1)  ||  tpff  lla.iT,  +  II  <Pn  lli.fi) 

since  ci(pn,iP)  =  0  for  F  =  0  and  supp  ui  C  B.  Recalling  (4.2)  with  t  —  p n  and  its 
corollary  (4.2'), 


II  n(upN)  ||2.ir,<  C{p  m‘n{U(7  1)}  ||  <pN  (J,  Ki  +  ||  <pN  ||_fc|Kl). 

Adding  the  two  upper  bounds  yields  (4.3). 


□ 
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Lemma  4.4  Under  the  hypothesis  of  the  previous  lemma,  with  0  <  k  <  5 

II  <Pn  IU./c,, <  C  ||  <pN  ll-jt.if,  (4.4) 

with  C  depending  on  K0,Ki,k  and  pt(<p)  =  7  >  2  —  y. 

Proof:  Let  K0  CC  Kx/i  CC  •••  CC  CC  K\  be  unions  of  elements;  i  was 

defined  just  below  (4.2).  Lemma  4.3  applies  to  each  pair  <  j  <  k.  So 

for  7  <  3/2, 

II  <PN  IU, C(p~2^~1')  II  <pN  ||2,K-1)y+1)  +  II  <PN  II /f,(y+l)) 

Iterating  from  j  =  0  to  A:  +  1,  we  get 

II  <PN  II 2,JFCo  —  C(p-2^~1Hfc+1)  ||  <pN  ||2iK-.(t+I)  +  ||  <pN  || 

Let  if,(fc+i)  CC  K'  CC  Ki.  But,  by  Schmidt’s  inequality  [10], 

p-2(fc-2)  ||  if>N  ||2  K..(t+i)  <  p"2(fc-2)  ||  <pN  §2  K, 

<  C  ||  <PN  II ~k,K' 

<  C-  II  <pN  H-fc./f,, 

since  (A:  +  1)(7  —  1)  >  k  —  2  iff  7  >  2  —  *|y  and  this  last  expression  <  3/2  iff  k  <  5. 
The  case  7  >  3/2  follows  in  a  similar  way,  again  we  need  k  <  5. 

□ 

Next  the  local  version  of  the  main  result  of  the  section. 


Lemma  4.5  Let  Kq  CC  K\  be  unions  of  elements,  <p  6  if‘(.Ki),^/v  £  Yn  and  k  G  N. 
Then  for  2  <  s  <  1  +  7, 

II  e  |U„<  C(p-*-+»  +  ||  e  IUk.)  (4.1) 

Proof:  Let  Kq  CC  K1/3  CC  K2/3  CC  K\  and  w  =  1  on  Ki/$,w  G  Co^ifys).  Then 
||  w<p  -  1 Ua.if,  <  C  inf  ||  u<P  -  IU.it, 

S>£Yn{K  j) 

<  Cp-2^ 


and 


\<P~<Pn  1 1 2 ,  iCo  —  1 1  IP  ~  *r(u/£>}  IU.k-o  +  II  "'(wip)  “  *Pn  tU.iTo 


14 


where  ci(ipN  —  n{uj(p),ip)  =  0,  Vip  £  Yn(K x/3)  so  that  lemma  4.4  applies  with  Kx/ 3 
replacing  Kx , 


||  <pN  -  ir(w<p)  ||a>*0  <  C7  ||  —  7r(urv>)  ||_fc,#cl/s 

<  C{\\  e  ||_fc.*I/3  +  \\u<p-  n(u<p)  H-*.*,,,) 

<  <7(11  e  11.*^,  +  ||  wip  —  7r(w«p)  H-*.*,) 

and  collecting,  we  get 


II  e  ||2,*0  -  ||  <P~  <Pn  lk*0 

<  0(11  uip-  n[u<p)  H-fc.K,  +  ||  e  (I-*.*,) 

<  C(p~2^+  ||  e  ||_*.*,) 


which  is  (4.5)  when  s  =  2.  Let  s  >  2,Ko  CC  ifi/3  CC  K2/ 3  CC  Kx  as  before.  Then 
for  some  4>  €  Yn  , 


II  «  II., K.  <  II  V  -  <P  ||..K,  +  II  <l>  -  <pn  ||., K,,. 

<  II  P  -  i>  If,  +P!("”  II  *  -  f  llw,  +P!,"’>  II  e  ||!lKl 

<  C(p-,11"+1 1  +  p=(— =)-*(->—+»  +  pX-«)  ||  e  ||.,iK.) 

<  c(p-!<’-'+»+p1'-,l||e||.,,K1) 

which  is  (4.5). 

The  main  result  generalizes  (4.5)  to  hold  for  flo  CC  fli. 


□ 


Theorem  4.6  Let  fl0  CC  fix  CC  R2,  <p  €  tf‘(Hi),  <pN  e  YN  and  k  £  N. 

2  <  s  <  1  +  7, 


II  *  ll.,n„<  C(p-’<’-«>  +  p*<-»  ||  e  II-,, n,)- 


Then  for 
(4.6) 


We  may  think  of  fix  =  fi  or  fli  CC  fi  depending  on  the  regularity  of  <p. 

Proof:  Cover  fi0  with  a  finite  number  of  ff0(xj),  i  =  centered  at  x,  £  fl0 

with  sidelengths  =  |dist(flo,dni).  Let  Kx(xi)  be  the  same  except  with  sidelengths 
=  dist(fi^,dfii).  For  s  >  0,  we  have  from  the  previous  lemma  that 

II  f  ||.,k.(„)  <  c,(p-!<’-*"  +  p’<-»  ||  f  H-w.w) 

<  Ci(p-,<’"’+1>  +  p’(— *>  ||  e  ||-i,n.) 


which  yields  (4.6). 


□ 


Corollary  4.7  With  the  same  hypotheses  as  in  Thm.  4.6, 


II  *  IU<  (4.7) 

Proof:  For  k  >  0,  2(s  —  2)  +  k  —  27  >  —2(7  —  5  +  1)  and  applying  cor.  3.6  yields 
(4.7). 

□ 

To  obtain  (4.7)  it  is  thus  sufficient  to  use  k  =  0,  i.e.  no  restrictions  are  needed  on 
7  other  than  7  >  0,  see  lemma  4.4.  This  provides  quasi-optimai  error  estimates  for 
s  >  2. 

Remark  4.8  For  F  ^  0,  we  can  modify  the  lemmas  to  reflect  <px  satisfying  instead 

E[xl})  d=  cx(e, t/»)  =  Ci(<p~(pN,rp),  V0  G  YN{KX)  (4.8) 

where  E  G  \H^(KX)}\  For  5  >  -2,  s  G  Z,  let  ||  E  ||_,>n=  supu6W,+2(n)  \E(rp)\/ 

II  v  ||,+2,n-  If  <Pn  satisfies  (4.8),  lemma  4.1  should  state  instead, 

||e||-«,jf0  <  C{p-<\\e\\itKl  +  ||«||— •— i.itj  +  P-4||^||2.k,  +  ||.E||_,,k,); 

and  lemma  4.2, 


||e||o.Kn  <  C(p  4||el|2,if,  +  llcll-fe,iC,  +P  4||-^l|2.ir,  +  ll^llo.A-,)- 

The  remaining  lemmas  and  theorems  hold  unchanged.  In  the  proofs  one  has  to  carry 
along  and  estimate  terms  ||J5||r- 

□ 

Remark  4.9  It  is  not  always  necessary  to  have  H0  CC  Hi  or  K0  CC  Ku  since  it  is 
at  times  possible  to  extend  about  parts  of  the  boundary  by  the  methods  mentioned 
in  remarks  3.2  and  3.3  (or  the  smoothness  of  the  solution  otherwise  warrants  this). 
In  such  cases,  we  may  use  the  estimates  up  to  that  part  of  the  boundary. 

□ 


5.  A  projection  method. 

We  now  discuss  how  to  recover  the  hydrostatic  pressure  utilizing  that  it  satisfies  the 
Poisson  problem  (1.6).  We  shall  see  that  there  seems  to  be  a  fundamental  interplay 
between  the  lacks  of  stability  and  regularity. 

First  suppose  P  is  smooth,  i.e.  P  G  H*  n  W  for  some  s  >  1.  We  will  as  before 
assume  that  the  regularity  of  P  is  dictated  by  that  near  the  boundary  of  0.  We  will 
then  use  the  following  weak  formulation 
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Find  PeH'nW  such  that  Vg  €  H *(n)  n  W , 

*(P,«)=f(U,«)  +  (F,V*)  (5.1) 

where  the  forms  in  (5.1)  are  defined  by 

di(P,q)  =  [  VP  •  Vg  dx, 

•'n 

jf(U,  g)  =  —  /  AU  •  n  9  ds.  (5.2) 

■'an 

In  order  to  discretize,  we  use  instead  <7(U*,g),  which  is,  of  course,  why  we  needed 
higher  order  error  estimates  for  ||U  —  TJ tv || ,  since  for  s  >  0, 

|| AU  •  n  -  AU*  •  n|U_3/2.an  <  C-||U  -  U*||1+g,n.  (5.3) 

Note  that  it  is  possible  (and  sometimes  desirable)  to  pose  (5.1)  locally,  i.e.  over 
a  union  of  elements,  rather  than  over  all  of  Q.  We  may  discretize  selecting  WN  — 
Pp_1  n  W .  We  will  estimate  in  turn  each  of  the  two  terms  on  the  right-hand-side  of 

||P  -  Pn\\i,Ko  <  ll-P  -  -Puwlli.x,,  +  ||Puw  -  -fWlli./r,, 

where  P\jN  satisfies  di(PuN,g)  =  </(U*,g),  V<?  6  W*.  Let  q  =  P  -  P\jN ,  then  since 
dj  is  symmetric  and  coerces  the  1-norm, 

\\P~PvJla  <  ||P-i>uJ|8/2-.,an||A(U-U*).n||._3/2>an 
<  C||P-PuJ|2_<in||U-U*||1+,(n 

which  with  s  =  1  and  (4.7)  gives 

\\P  -  Fuji.  < 

For  the  second  addend,  P  inherits  a  singularity  from  U  in  (1.6),  P  =  r',-20(0)  so 
that  at  worst, 

l|J=W  -  J=o„lli  < 

Hence, 

WP-PAiSCp-^-'K  (5.4) 

Example  5.1  Consider  the  “driven  cavity”  flow  problem 

-AU  +  VP  =  0  in  J2  =  (-1,1)* 

V  •  U  =  0  in  R 

U(x,  1)  =  (1  -  x*)fl  ^  j  for  -  1  <  x  <  1 
U  =  0  elsewhere  on  dR. 
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where  0  >  0.  According  to  Thm.  6.2  of  [2|,  the  last  three  equations  has  a  solution  in 
\Hl+0~t{R)\2  and  no  solution  is  in  [Hl+^(R)\2 .  This  then  carries  over  to  hold  for  the 
solution  of  the  full  system  by  known  regularity  results  for  Stokes  problem  [18],  [12]. 

□ 


Proposition  5.2  Let  f!  =  (—1,1)*.  WN  =  Pp_x  n  W  and 

Find  PN  G  Wn  such  that  Vq  G  Wn ,  /_  _v 

d1(PN,q)=g(UN,q)  1  j 

where  dt,g  and  U n  G  (#o)2  have  been  defined  above.  For  the  driven  cavity,  (3  >  l, 
we  obtain 


||P-P/v|kn  <  Cp-^-V  and 

||  P  —  P/v||o,n  <  Cp"2^-1)  (5.6) 

Note  that  the  second  half  of  (5.6)  corresponds  exactly  to  the  rate  of  approximation 
of  the  velocity  in  H1. 

Proof:  For  the  square,  the  odd  extension  procedure  works  for  U,U^  G  {Hq)2  and 
U|(_U)X{1}  G  PTq ^ ( ( —  1 » 1 )  x  {!})  in  such  a  way  that  the  smoothness  properties  are 
preserved  and  the  norms  of  the  extended  functions  remain  the  same  as  those  of 
the  original  functions.  Applying  (5.4)  yields  the  1-norm  estimate  in  (5.6).  To  get  the 
second  half,  we  employ  Aubin-Nitsche’s  construction:  Let  6  =  P  —  Pn  and  x  £  HoW 
solve  the  auxilliary  problem 


di(x,^)  =  M),  W >ew. 

Then  Hxlh  <  Cpllo  because  fl  is  convex.  ||5||g  =  (t,S)  =  d i(b,x)  =  di(S,x  -  t/>)  < 
C||tf||i||x  -  V'lli  <  Cp~2||$||i||£||0  as  we  can  choose  0  so  that  ||x  -  V'lli  <  C'p~2||x||j. 

□ 

Remark  5.3  Thus  we  are  able  to  explain  why  it  should  be  possible  to  get  optimal 
convergence  rate  also  for  the  pressure  as  observed  in  |16|  for  0  =  2.  Pn  defined 
in  (5.5)  solves  the  same  problem  as  in  the  method  from  [16]  since  we  observe  from 
Prop.  2.1  that  3p0  G  Z+,  such  that  for  p  >  p0,  <>(U^,g)  =  0,  Vg  G  WN  and 
6(v,P/v)  =0,  Vv  G  Zn,  rotations  of  streamfunctions.  Now  let  v  =  Vg  in  (1.3)  and 
integrate  by  parts  to  get  (5.5). 

□ 

Remark  5.4  With  P  smooth,  it  is  also  possible  to  introduce  mixed  methods  to  solve 
for  Pn  as  in  [26]  and  use  either  Raviart-Thomas  or  Brezzi-Douglas-Marini  elements 
to  get  almost  optimal  order  convergence  for  VP. 
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□ 

Remark  5.5  We  can  interpret  (5.1)  as  local  averaging  bringing  this  method  into  the 
same  general  frame  of  ideas  as  that  of  [7]  or  [17]. 

□ 

For  a  pressure  P  which  is  not  assumed  smooth  globally,  P  e  Ha  fl  W  for  some 
s  >  0,  the  above  procedure  can  be  employed  locally  on  an  element  or  a  patch  of 
elements  in  the  interior  of  fl,  where  P  could  be  assumed  smooth.  This  would  yield 
an  interior  optimal  rate  approximation  to  P. 

For  P  still  not  smooth  one  might  be  tempted  to  modify  (5.1)  to  some  other  weak 
form  of  (1.6)  presupposing  minimal  smoothness  on  P. 

Example  5.6  In  stead  of  (5.1)  take  the  following  problem  on  fl  =  ( — 1,  l)2: 

Find  P  £W  such  that  Vg  6  H2  n  {q  :  §*  =  0  on  50}  ,  . 

d2(P,q)=g(U,q)  1  j 

where  the  bilinear  form  d2  is  defined  by 

d2 (PiQ)  =  ~  j  P&qdx  (5.8) 

■'n 

Note  that  for  the  discrete  form  of  (5.7),  P ^  would  be  identical  to  that  of  (5.5)  since 
0  =  •/an^vf*  ds  =  d2(PN,q)  +  di(PN,q)  for  q  6  WN.  So  one  might  expect  to  get 
optimal  rate  0-norm  estimates  for  P  —  Pjq,  since  also 


inf 

Pew 


sup 

qeH*,  |j|=o  on  an 


MP,q)  >c^n 
Mollill*  - 


because  Vp  €  W  \{0},3g*  €  2/2n{j£  =  0  on  50}  such  that  A q*  =  p,  ||g*||j  <  C||p||0- 
Nevertheless  one  will  fail  to  get  optimal  rate  of  convergence  since  d2  is  p-unstable, 
see  (2.2.28)  in  [3]  for  one  dimension. 

□ 


6.  A  method  using  potentials 

We  will  demonstrate  pressure  recovery,  still  using  the  Poisson  equation  (1.6),  but  now 
in  a  way  so  that  we  get  the  hydrostatic  pressure  evaluated  at  an  interior  point  of  0. 
Let  F  =  0  so  that  P  satisfies,  for  any  fii  C  0  with  50i  sufficiently  regular, 

A P  =  0  in  Oi 

=  ( AU)  •  n  on  50j 


(6.1) 
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Then  ,  if  dCli  G  C3, 

F(x)  =  -\ak\U0S  |X-^|(AU'”H>,)  <<•(*>  <62> 

and  PuN  is  obtained  by  replacing  U  in  (6.2)  by  U^. 

Corollary  6.1  If  P\ jN  is  evaluated  exactly  on  dCli,  smooth,  x  G  flo  CC  fi|  CC  fi 
according  to  (6.2)  with  U#  in  stead  of  U,  and  U  G  Hl+t(U1)  for  some  s  >  0  and  of 
singularity  strength  7  —  1,  then  for  k  G  N, 

| P  ~  AjJ(x)  <  C{ih)p-k~2^  (6.3) 

Note  this  represents  a  superconvergence  result  for  the  p  version  in  the  interior  of 
elements. 

Proof:  From  (5.3),  (3.3),  and  the  fact  that  the  single  layer  potential  is  smooth  away 
from  <917  j. 

□ 

Remark  6.2  It  might  be  possible  to  deal  with  dfli  — *■  <917  by  using  explicit  knowledge 
about  U  near  corners  of  17.  Such  explicit  knowledge  can  be  generated  a  priori,  e.g 
using  the  methods  of  P.  Papadakis,  [21]. 

□ 

Remark  6.3  It  might  be  desirable  to  evaluate  the  right  hand  side  of  (6.2)  with 
replacing  U  by  quadrature. 

□ 

Remark  6.4  We  note  that  this  procedure  can  be  performed  without  smoothness 
assumptions  other  that  U  G  Hl+,{ Hi),  i.e.  for  P  G  H* ,  s  >  0. 

□ 

Other  such  extraction  procedures  have  been  investigated  for  nearly  incompressible 
plane-strain  elasticity  in  (27).  The  two  continuous  problems  are  connected  as  shown 
in  [29],  Chap.  6. 


7.  Concluding  remarks 

We  have  introduced  two  methods  for  recovering  the  hydrostatic  pressure  in  Stokes 
problem.  The  necessity  for  this  was  demonstrated  in  [16].  We  derived  negative  norm 
and  interior  estimates  enabling  us  to  prove  optimal  convergence  rate  in  the  interior 
elements  and  for  special  geometries  and  boundary  conditions  up  to  the  boundary. 
Thus  we  were  able  to  explain  the  ability  of  the  mixed  method  introduced  in  [16]  to 
approximate  P  at  optimal  rates  when  the  boundary  data  were  sufficiently  smooth. 


20 


References 

[1]  R.  A.  Adams.  Sobolev  Spaces.  Academic  Press,  1975. 

[2]  D.  Arnold,  L.  R.  Scott,  and  M.  Vogelius.  Regular  inversion  of  the  divergence 
operator  with  Dirichlet  conditions  on  a  polygon.  Arm.  Sc.  Norm.  Super.  Pisa, 
Cl.  Set.,  IV  Ser.,  15(2):169-192,  1988. 

[3]  D.  N.  Arnold,  I.  BabuSka,  and  J.  Osborn.  Finite  element  methods:  principles 
for  their  selection.  Comput.  Methods  Appl.  Mech.  Engrg.,  45:57-96,  1984. 

[4]  I.  BabuSka.  Error  bounds  for  the  finite  element  method.  Numer.  Math.,  16:322- 
333,  1971. 

[5]  I.  BabuSka  and  M.  Suri.  On  the  optimal  convergence  rate  of  the  p-version  of  the 
finite  element  method.  SIAM  J.  Numer.  Anal.,  24:750-776,  1987. 

[6]  I.  BabuSka,  B.  Szabo,  and  I.  N.  Katz.  The  p  version  of  the  finite  element  method. 
SIAM  J.  Numer.  Anal.,  18:515-545,  1981. 

[7]  J.  H.  Bramble  and  A.  H.  Schatz.  Higher  order  local  accuracy  by  averaging  in 
the  finite  element  method.  Math.  Comput.,  31(137):94-111,  1977. 

[8]  F.  Brezzi.  On  the  existence,  uniqueness  and  approximation  of  saddlepoint  prob¬ 
lems  arising  from  Lagrangian  multipliers.  RAIRO,  8:129-151,  1974. 

[9]  M.  Dauge.  Stationary  Stokes  and  Navier-Stokes  systems  on  two-  or  three- 
dimensional  domains  with  corners,  part  I:  Linearized  equations.  SIAM  J.  Math. 
Anal.,  20(l):74-97,  1989. 

[10]  M.  Dorr.  The  approximation  theory  for  the  p-version  of  the  finite  element 
method.  SIAM  J.  Numer.  Anal.,  2 1  (6) :  1 1 80— 1 207 ,  1984. 

[11]  D.  Gottlieb  and  S.  Orszag.  Numerical  Analysis  of  Spectral  Methods:  Theory  and 
Applications.  Volume  26  of  CBMS-NSF  Regional  Conf.  Series  in  Applied  Math., 
SIAM,  1977. 

[12]  P.  Grisvard.  Elliptic  Problems  in  Nonsmooth  Domains.  Pitman,  1985. 

[13]  T.  J.  R.  Hughes,  L.  P.  Franca,  and  M.  Balestra.  A  new  finite  element  formulation 
for  computational  fluid  dynamics:  V.  Circumventing  the  BabuskarBrezzi  condi¬ 
tion:  A  stable  Petrov-Galerkin  formulation  of  the  Stokes  problem  accomodating 
equal-order  interpolation.  Comput.  Methods  Appl.  Mech.  Eng.,  59:85-99,  1986. 


21 


[14]  J.  Douglas,  Jr.  and  J.  Wang.  An  absolutely  stabilized  finite  element  method  for 
the  Stokes  problem.  Math.  Comput.,  52(186) :495-508,  1989. 

[15]  S.  Jensen.  An  H interpolation  result.  Submitted  to  SIAM  J.  Math.  Anal. 

[16]  S.  Jensen  and  M.  Vogelius.  Divergence  stability  in  connection  with  the  p  version 
of  the  finite  element  method.  RAIRO,  Modelisation  Math.  Anal.  Numer.,  1990. 
To  appear. 

[17]  C.  Johnson  and  J.  Pitkaranta.  Analysis  of  some  mixed  finite  element  methods 
related  to  reduced  integration.  Math.  Comput.,  38(158) :357-400,  1982. 

[18]  R.  B.  Kellogg  and  J.  E.  Osborn.  A  regularity  result  for  the  Stokes  problem  in  a 
convex  polygon.  J.  Funct.  Anal.,  21(4):397-431,  1976. 

[19]  J.  A.  Nitsche  and  A.  H.  Schatz.  Interior  estimates  for  Ritz-Galerkin  methods. 
Math.  Comput.,  28(128):937-958,  1974. 

[20]  J.  Osborn.  Regularity  of  solutions  of  the  Stokes  problem  in  a  polygonal  domain. 
In  B.  Hubbard,  editor,  Symposium  on  Numerical  Solutions  of  Partial  Differential 
Equations  III,  pages  393-411,  1975. 

[21]  P.  Papadakis.  Computational  aspects  of  the  determination  of  the  stress  intensity 
factors  for  two-dimensional  elasticity.  PhD  thesis,  Univ.  of  Maryland,  1988. 

[22]  E.  M.  R0nquist.  Optimal  spectral  element  method  for  unsteady  three-dimensional 
incompressible  Navier-Stokes  equations.  PhD  thesis,  MIT,  1988. 

[23]  L.  R.  Scott  and  M.  Vogelius.  Conforming  finite  element  methods  for  incompress¬ 
ible  and  nearly  incompressible  continua.  In  Lectures  in  Applied  Mathematics, 
pages  221-244,  AMS,  1985. 

[24]  S.  Serbin.  A  Computational  Investigation  of  Least  Squares  and  other  Projection 
Methods  for  the  Approximate  Solution  of  Boundary  Value  Problems.  PhD  thesis, 
Cornell  Univ.,  1971. 

[25]  M.  Suri.  The  hp  version  of  the  finite  element  method  for  elliptic  problems  of 
order  2  l.  RAIRO,  Modelisation  Math.  Anal.  Numer.,  24:107-146,  1990. 

[26]  M.  Suri.  On  the  stability  and  convergence  of  higher  order  mixed  finite  element 
methods  for  second  order  elliptic  problems.  Math.  Comput .,  54(189):1-19,  1990. 


*  “ 


22 

[27]  B.  Szabo,  I.  Babu§ka,  and  B.  K.  Chayapathy.  Stress  computations  for  nearly 
incompressible  materials  by  the  p-version  of  the  finite  element  method.  Int.  J. 
Numer.  Methods  Eng.,  28(9):2175-2190,  1989. 

[28}  B.  A.  Szabo.  PROBE:  Theoretical  Manual.  Noetic  Technologies  Corp.,  St.  Louis, 
Missouri,  1985. 

[29]  R.  Temam.  Navier-Stokes  equations.  North-Holland,  1977. 

[30]  M.  Vogelius.  A  right-inverse  for  the  divergence  operator  in  spaces  of  piecewise 
polynomials  -  Application  the  the  p  version  of  the  finite  element  method.  Numer. 
Math.,  41:19-37,  1983. 


